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We use a path integral approach for solving the stochastic equations underlying the financial 
markets, and we show the equivalence between the path integral and the usual SDE and PDE meth- 
ods. We analyze both the one-dimensional and the multi-dimensional cases, with point dependent 
drift and volatility, and describe a covariant formulation which allows general changes of variables. 
Finally we apply the method to some economic models with analytical solutions. In particular, we 
evaluate the expectation value of functionals which correspond to quantities of financial interest. 



INTRODUCTION 



The starting point of our analysis is that quantities as the stock prices, the option prices, and the interest rates 
satisfy differential stochastic equations (SDEs), i.e. ordinary differential equations with a superposed white noise. 
Such equations are called Langevin equations, and they are extensively used in the financial literature. The solutions 
of such equations are usually obtained by solving the associated partial differential equations (PDEs). 

In this paper we want to describe an alternative approach based on the path integral formulation. The notion of 
path integrals, also called Wiener integrals in stochastic calculus, and Feynman integrals in quantum mechanics, is 
known from a long time. A proper mathematical definition of the Wiener integral can be found in the original works 
of Wiener jjj , and Kac , while the quantum mechanical analogous has been introduced by Feynman || . 

The importance of this formalism lies in the possibility of employing powerful analytical and numerical techniques, 
developed in physics, for solving the usual problems of option pricing. Some attempts of using this approach in finance 
have been described in recent literature (see, for instance, Ref. Q). Here we discuss the path integral formulation in 
a general manner, and, as examples, we solve some well known economical models. Furthermore, in a forthcoming 
paper, we will describe some numerical methods. 

In section |, we show the general one-dimensional formalism. In section [n|, we extend the formulation to the 



multi-dimensional case. In section III, we discuss a covariant formulation which is necessary to perform a general 



transformation of variables. In section IV, we define the expectation value of a general functional. Finally, in section 
[y|, we give some analytical results. 

I. TRANSITION PROBABILITY AND PATH INTEGRAL FORMALISM IN ONE DIMENSION 



A. Langevin equation and discretization problem 

In 1908 the french physicist Paul Langevin |5| wrote down a differential equation containing a Gaussian white noise 
coefficient, /(t), 

x(t) = a(x,r) + a(x,T)f(r). (1) 

where x(t) is a stochastic process to be determined. Such white-noise-driven differential equations are often used 
in physics and chemistry and are the oldest form of SDEs. The Langevin equation written above does not define 
univocally a stochastic process, and it has to be supplemented with an additional interpretation rule (see, for instance, 
Refs. [^-||). This is related to the ambiguity in the discretizations of this equation. We may remove this ambiguity 
by introducing additional structure into the Eq. (|l]), and discretizing the Langevin equation as 

Ax = a(y + (Ax, t)At + a(y + rjAx, t)Aw, (2) 
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where t is the initial time, At the time step, and y — x(t). If we expand in Ax and recall that Aw ~ 0{\fAt), we 
find, Ax = o(y, t) Aw + O(At), and, to the leading order in At, 

do 

Ax = a(y, t)At + v o{y, t) — (y, t) Aw 2 + a(y, t) Aw. (3) 

Now, since in the framework of stochastic calculus the following equality holds, 

Aw 2 =At, (4) 
(where we used the usual symbol |10|j , =), the Eq. (|J) becomes 

Ax = A(y,t) At + o(y,t) Aw, (5) 

where 

A(x,t) = a(x,r) + i] <t(x,t) — (x,t). (6) 

The Langevin equation, written in this form, describes a well defined stochastic process, and different r/'s correspond 
to distinct processes. For example, r\ = and r\ = — correspond to the Ito and the Stratonovich interpretations, 
respectively. Note that the Eq. (|5|) does not depend on £. In conclusion, 

I. A stochastic differential equation is well defined only if both a continuous expression and a discretization rule 
are given. 

From now on, we will always write the underlying stochastic equation in the discretized form (^) , understanding that 
the continuous limit must be taken. Therefore the stochastic process will be defined when the functions A(x, r) and 
o(x,t) are given. 

Unfortunately there is still an ambiguity. In fact, even if a stochastic process is well defined, the SDE describing 
such a process is not univocal. Let us consider the following family of equations, depending on the parameter, n, 

do 

x{t) = a(x,r) ~ k o(x,t) — [x, t) + o(x, r)/(r). (7) 
By expanding in Ax, as above, we obtain 

do 

Ax = [a(y, t) + (r)- k) o(y, t) -£(y, t)} At + o(y, t) Aw. (8) 

Therefore, if we choose the appropriate discretization rule for each parameter, k, i.e. by fixing r\ such that r\ — k is 
fixed, we can describe the same process by different SDEs. In other words, 

II. Many different continuous expressions for the SDE, with the appropriate discretization rule, define the same 
stochastic process. 

We want to stress that all economic and financial applications of stochastic calculus have used so far the Ito 
prescription because it gives the coefficients of the SDE a simpler and meaningful interpretation (in particular, the 
drift coefficient, a(x, r), appears directly). On the other hand, the Stratonovich prescription, for example, allows 
to employ the usual rules of calculus instead of the more complex Ito calculus. However we have seen above that 
any SDE with the Ito prescription is equivalent to some other SDE with the appropriate prescription. This more 
general approach is needed in order to describe the connection with the equivalent difficulties which arise in the path 
integral formulation, where we cannot avoid a general treatment to explain the analytical and numerical techniques 
of computation. 

Finally, we observe that an alternative description of the stochastic process is given by the following equivalent 
partial differential equations, which can be derived in an unambiguous way from the Langevin equation (||) (see, for 
example, Ref. 0]): 



Kolmogorov's backward equation 



^p( x ,T\y,t) = ^A{y,t)^-+ l -o 2 (y,t)^\ ,,[.,: T | ,/. /): (!)) 
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• Kolmogorov's forward equation or Fokker-Planck equation 

-^p(x,T\y,t) = j-A A( a:> T) + ~(7 2 (a; > T)| p(x,T|y,i). (10) 

The solution of these equations is the transition probability function, p(x,T \ y,t). Note that the ambiguities of the 
continuous Langevin equation (Q) are not present here for the following reasons: the first ambiguity is solved once 
the functions A(x,t) and cr(x, t) are given; the second one, which is connected to the differential operator ordering, 
is fixed once this ordering has been fixed. 



B. Stochastic differential equation and short-time transition probability 

The Brownian motion can be seen as the convolution of an infinite sequence of infinitesimal (short-time) steps. This 
constitutes the bridge between local equations and an integral formulation of the problem. Let us write the simple 
stochastic equation 

Ax = aAw. (11) 

If w is a Wiener process this equation define a Markov process with zero average and variance equal to a 2 At. The 
corresponding short-time transition probability is given by 



p( X , t + At | y, t) = y exp { JLJL } . (12) 

In general, we do not have an explicit expression for the short-time transition probability corresponding to the 
stochastic equation (^|). However we can write the following general expression, which is correct up to 0(At) 1 



P (x, t + At | / 1 , exp ~ {x - V - - % M) , (13) 

V 2tt At a(y,tf P \ 2 At a(y,t) 2 J 

where A(x, r) is given in Eq. (0). A proof of this is reported for completeness in Appendix |a|. The equation ( |l3| ) gives 
a prescription to write the solution of a Fokker-Planck equation in the form of a convolution product of short-time 
transition probability functions. 



C. Finite time transition probability and path integral 



The finite time transition probability can be written as the convolution of short-time transition probabilities 

/+oo r+oo 
. . . dxjydxiy^i . . . dx\ p(x, T | xn, T — At) 
-oo J — oo 

p(x N , T - At | x N -i,T - 2At) . . . p( Xl ,t + At | y, t), 



with At — 



T-t 
N 

p(x,T | y,t) 



Therefore, by substituting the expression (|13|) in the previous equation, we get 



+oo r +oo N N+l , .J 

i=l 3=1 



oo *J — oo 



y / 2nAt 



with 



So(xj, tj] Xj-i, tj—x) 



2 a(xj-i,tj-i) 



(Xj Xj — l) 

At 



A(xj-i,tj-i) 



At 



(14) 



(15) 



(16) 



(the meaning of the subscript, 0, will be explained in the next section). We can interpret — 3 ^ as a mean 

velocity in the time interval At, and cr(xj-i, tj—i) as a mass. Then a Lagrangian structure appears explicitly. 



Hence, the limit for N — > oo of the RHS of Eq. (|15|) can be formally written as 
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p{x,T | y,t) 



x(T)=x - \ f T 

V[<j{x,t) 1 x(t)\ cxp<- / L [x{t),x(t);t] dr 

x(t)=y I Jt 



(17) 



The functional measure, means summation on all possible paths starting from x(t) — y and arriving 

at x(T) — x. The integral to the exponent must be interpreted as the limit of the discrete summation for N — > oo, 
and the function Lq[x(t), x(t); t] is a Lagrangian function defined by 



Lo\x,x:t] = 



1 



2 <t{x,t) 2 



[x — A(x, t) 



(18) 



The RHS of Eq. @ is called "path integral" . Since in the limit N — * oo only O(At) contribute to the integral, the 
formal expression in Eq. (|l^) gives the exact finite time transition probability. In conclusion, there is a complete 
equivalence between the differential stochastic equation, the Fokker-Planck equation and the path integral approach. 

Note that in the limit N — » oo distinct expressions for the short-time transition probability, equal up to O(At), give 
rise to the same functional expression for the path integral; therefore orders greater than O(At) do not have any effects 
from an analytical point of view. However, if the path integral is solved numerically, a more accurate approximation 
for the short-time transition probability can be useful. A general procedure to obtain an approximation to any order 
is described in Appendix [E|. 



D. Path integral and discretization problem 



We note that the path integrals present exactly the same ambiguities of the stochastic differential equations. In 
particular, it is clear from its definition that the formal expression (|l7j ) is well defined only if a discretization rule, 
i.e. a short-time transition probability, is also given, since in general different discretization rules give rise to different 
results. Therefore we can say, 

III. A path integral is well defined only if both a continuous formal expression and a discretization rule are given (cf. 
with I). 

A further ambiguity is due to the fact that different Lagrangians, discretized by the appropriate rule, give rise to 
equivalent expressions for the short-time transition probability. For example let us consider the Lagrangian 



L 1/2 [x,x;t] = L [x,x;t] + - d x A(x,r), 



with the discretization rule given by 
p{x,T \y,t): 



oo «/ — OO 



N N+1 a ( x - t ) _1 

Y[dxi JJ 3 --j== exp{-Si/2(xj,tj;Xj-i,tj-i)} , 



i=i j=i 



V2nAi 



where 



Sl/2(Xj, tj] Xj-i, tj-i) — 



At 



2 a{x j - 1 ,t j - 1 f 



(xj Xj-i) A/ Xj + Xj—i tj + tj_x ^ 
At ( 2 ' 2 [ 

By expanding in series of Ace and At up to O(Ai), we obtain 



^1 d x A( Xj+Xj ^ 1 tj + t ^ x 



Sl/2( x j> tj] Xj^i, t j — 1 ) 



At 



2 fr^j-i.tj-i) 



{Xj Xj—x) 

At 



At Ax 
2 I a[Xj-i,tj-i) 



— A[Xj_i, tj_Xj 

2 d x A(xj-i,tj-i). 



Then by using the following identity, analogous to the (Q), 

Aa; 2 = o-(xj-i,tj-if At, 



(19) 



(20) 



(21) 



(22) 



(23) 
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the Eqs. (|2(]) and become equal to the Eqs. ( |l5| ) and ([l6]). The formal expression of the path integral is 

, f x(T)=x ( ,T } 

p(x,T\y,t)= Viaix^y^ir)} expl- L 1/2 [x(t),x(t);t] dr) . (24) 

J J x (t)=y { Jt J 

We will call pre-point the path integral formulation given by the Eqs. jl5|), (|l|), (17), and (|l8|), and mid-point that 
given by the Eqs. @, @, (|l|), and @. In conclusion, 

IV. Many different continuous formal expressions for the path integral, with the appropriate discretization rule, 
define the same stochastic process (cf. with II). 

The pre-point formulation is usually simpler to perform numerical computations, while the mid-point one has some 
advantages to carry on analytical calculations, as we will see in the following. From now on, if not explicitly specified, 
we will use the mid-point formulation, but for the sake of simplicity we will omit the subscript 1/2. 



Let us consider the case 



E. Gauge transformation 



L[x, x; r] = -^-^ [x - a{x)f + i d x a{x), 



(25) 



where we have taken a(x,r) = a = constant, and A(x,t) — a{x). The finite time transition probability is given by 



p{x,T | y,t) 



:(T)=x 



V[a 1 x(r)]exp< 



dr 



(26) 



This expression contains a coupling between the derivative of the stochastic variable, x, and the term a(x). Let us 
perform a gauge transformation by introducing a function 4>{x) such that 



a{x) = Tx< 



(27) 



then we have 



a{x) x dr 



dx 



dx = 



[<Kx)-<Ky)], 



(28) 



which is the same for all paths and depends only on the end points. 

It is important to note that the integrals in ( p8| ) are stochastic integrals (x is not defined for a Brownian path), 
and the result ( p8f ) is true only if we use the mid-point formulation (Stratonovich prescription). In fact, if we use the 

pre-point formulation, the term — d x a(x) in Eq. ( p6| ) should be dropped, while 



— ^ / a{x) x dr — I — dx 
v Jt ° Jt dx 



i [<P(x) - 0(y)] -I J ^<K<r)) dr 



(29) 



(Ito prescription). Of course the final result is the same in both cases. 

The term in Eq. (p8|) can be put out of the path integral, and Eq. becomes 



p(x,T | y,t) 



x(T)=x 



x(t)=y 



V[a 1 x(t)] exp ■ 



dr 



(30) 
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where 



Acj) = J a(x) dx. (31) 

In conclusion, the previous discussion shows that the convenience of adopting the mid-point formulation is related 
to the possibility of using the usual rules of integral calculus. 

II. MULTI-DIMENSIONAL CASE 

Until now we have discussed only the one-dimensional case. In the general case of a multi-dimensional Langevin 
equation in n dimensions, we have a straightforward generalization of it. Let us write the discretized Langevin 
equation in the following form, 

Ax" = A"(x(t), t) At + erf (x(i), t) Aw t (32) 
(the sum over repeated indices has been understood), where 

A"(x,r) =a"(x,r)+7 ? <(x,T) ^f(x,r). (33) 

The Eq. ( |32"|) corresponds to the following Fokker-Planck equation (in the sense that they describe the same stochastic 
process) 

4- p(x, T I y ,t) = A"(x,T) + — — — —7- — G^(x,T)l p(x,T I y,t), (34) 

where the differential operators in the right hand side act also on p, and 

G^(x,r) = cf(x,T)<(x,T). (35) 
The short-time transition probability is given by 



p(x J ,i i I xj-i.tj'-i) ~ (2^Ai)-"/ 2 ^G^-xttj-x) exp{-S* (x,,t J ;x,_ 1 ,^_ 1 )} (36) 

where 

5o(xj,tj-;Xj_i,tj_i) 

' G^fx^!,^!) [A^-A^x^t^Ai] [A^-A I/ (x i _ 1 ,t i _ 1 )At] , 



2 At 



(37) 



and x = y, x^r+i = x; while a path integral representation of the finite time transition probability is 

r x(T)=x _ f r T 

x(t)=y 

with 



p(x,T| y ,i) = // ( ) 2?[VGx(t)] expj- / L [x(t),x(t);t] drl, (38) 



L [x, x; r] = \ G MV (x, r) [i" - A M (x, r)] [x w - A" (x, r)] , (39) 

and the discretization rule given above. This expression for the transition probability corresponds to a well defined 
discretization procedure, i.e. the pre-point rule (cf. with Eqs. (15), (16), (jT^), and fll8|)). Another representation is 
given by the continuous formal expression 



p(x,T|y,t)= ff ( ] 2?[VGx(r)] expl- f L 1/2 [x(r),x(r);r] dr\ , (40) 

J J x(t)=y ! Jt I 
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with the Lagrangian 



1 d 

L1/2 [x, x; t] = L [x, x; r] + - — (x, r) , 
and the (mid-point) discretization rule 

p(xLj,tj | Xj-i,tj-i) ~ (27rAfT" /2 ^G(x~^t~i) exp{-S , 1/2 (x :) -,^;x i _i,t :) _i)} , 

where 



At 
2" Qo;'' 



Aa;1/ ^^ Xj +Xj_i ' 

"At ( 2 ' 2 j 



(41) 



(42) 



(43) 



(cf. with Eqs. @) and @ 



III. COVARIANT FORMULATION 



The path integral formulations given in the previous sections are not covariant. As a result, a general (non linear) 
transformation of variables in the path integral cannot be performed according to the usual rules of calculus. In the 
following we describe a covariant formulation of path integrals and we use this formulation to solve some 

specific problems. 

The tensor, G^ u , defined in Eq. (|35|), transforms under a general transformation of variables as a contravariant 
tensor (see, for instance, Ref. |l^]). If it is invertible - we will assume that this is always the case - it can be 
interpreted as the metric tensor of a Riemannian manifold. On the other hand, is not a contravariant vector, 
while the quantity 

1 q 

/i"(x, r) = i4"(x, t) ■■= — VG G^(x, r), (44) 



with 

G = G(x,r) = det G^(x,r), (45) 

transforms exactly as a contravariant vector. Therefore a covariant path integral representation of the solution of Eq. 
(E54|) is given by 



p(x,T|y,<)= ff * ) P[VGx(r)] expj- / £[x(r),x(r);r] dr \ , (46) 

J J x(t)=y Jt 



r x(T)=x _ f r T 

x(t)=Y 

where 

£[x, x; r] = \ [i" - h»] [x v - h"} + ^= ^ VG h» + i R, (47) 
The scalar R is the curvature 

R = G» v R\ Xv , (48) 

and 

_ d^s dT 
Al5ly dx u dx s ^ 1 A " 51 VTI 1 Sv ' 



dA _ IJ£ fw_ pr; r A _ pr; pA / 4q \ 

il A^f ^ x u5 x i/r? 1 uj/ 1 in ) V^-V 



r A _ 1 n sx ( 9G u s dG^s dG Vil , 
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As usual, the continuous formal expression ( fig ) must be interpreted as the limit of a discretized expression 

H-oo ,.+00 N N+l 

p(x, T\y,t)= ... J[ d n Xl J] A>(x,-,*j I x,-_i, tj-_i). 

•/— oo J — oo 8 - =1 



(51) 



3=1 



In general, however, even if a discretization rule is correct in a given coordinate system, i.e. the corresponding short- 
time transition probability satisfies the Eq. ( |34| ) up to O(At), after a transformation of variables it is not correct any 
more; but a covariant path integral formulation needs also a covariant discretization rule. A covariant discretization 
of the continuous expression (Ej) is given by |Il|jL2| ] 

1 



x [1 + G^WAx" + Mi—G^h" - Jj= JL VG h» - Lr) 
+(- (^G A „/i A ) - Lr^)^Ax v 
-\ (dpGvx) Ax»Ax"Ax x 

~{L (d„d„G a0 ) - L T ^ u r a x p )Ax^Ax"Ax a Ax' 3 /At 
+\{-G„ v WAx v + X - (dpGvx) Ax»Ax"Ax x /At) 2 ], 



(52) 



where tj = (tj + tj—\)/2 and all functions, if it is not explicitly specified, are evaluated at Xj-_i and at any time 
between tj—i and tj. 

With this definition of path integral, the formal measure in the expression (Eq) corresponds to the limit 



V[VG x(t)] = lira 



G(x,tjv+i) rr /G(xj,^) An 



At^o y (27rAi)" Al y (2ttA£) 



(53) 



Therefore, since the quantity \/G d n x is an invariant, after a change of variables the measure (53) transforms in the 
following way 



V[VG' x'(r)] = 



dx'„ 



dx v 



VWG x(r)], 



(54) 



where 



dx[ 



dx v 



is the Jacobian of the transformation calculated in the final point (note that the symbol, ', does not 



x.T 



mean derivation). 



A. The case of G M " = constant 



If G^ v is constant, the affine connections, T x , and the curvature, R, are equal to zero. Then the Lagrangian ( |47| ) 
becomes 



£[x, x; r] = \ [±" - A"(x, t)\ [x v - A"(x, r)] + \ JL A' 1 ( x. r > 



(55) 



(cf. with Eq. (]4l])), and the discretized expression (|52|), 

pfe,^ | x,-_i, *,-_!) ~ (2nAt)- n / 2 VG expj-^G^A^A^I 



x [1 + G^A^Ax" + At{--G^A»A V - - — A" 
+i Giv&fte'JL A x + i (-G^As") 2 ]. 



(56) 



8 



Note that the equation (|5^) is simply the expansion up to O(At) of the expression in (|42|). 

In conclusion, in the case of = constant, the covariant formulation becomes the mid-point one. Therefore, 
if we can make a change of variables whichgives a constant metric (flat space), the covariant path integral can be 
handled by the usual methods. In section |V| we will show that, if the resulting Lagrangian is a quadratic form, the 
path integral can be solved by using well known analytical methods. 



B. One-dimensional case 



In the one-dimensional case the covariant formulation simplifies significantly. The transformation rule for the metric 
tensor in one dimension is simply 



dx 



then G'n can be made equal to 1 everywhere by choosing 

x = J dx \J G\\. 



(57) 



(58) 



In this coordinate system the curvature R vanishes, as follows by Eqs. (|4S|) and (|50|); as a result, since R is invariant, 
it vanishes in every coordinate system. Then, by using the notation G 11 — ct(x,t) 2 , and G\\ = G = — ; — , the 



Lagrangian (f47| ) becomes 



where 



a(x, t) 2 



L[x, x; t] 



1 



2 a(x,r) 2 



[x - h(x,r)] 2 + 



ct(x,t) d h(x,r) 
2 dx a(x, t) ' 



h(x, t) = A(x, t) - - a(x, t) d x a(x, r). 



Finally, the transformation rule of the measure is 



Vla'ix'^y 1 x'(t)} = 



dx^_ 

dx 



Vlaix^r)' 1 x{t) 



(59) 



(60) 



(61) 



x,T 



C. Ito lemma 

Let us consider, for the sake of simplicity, the one-dimensional Langevin equation (JsJ) - The Ito lemma states that 
a function, z(x), of the stochastic variable follows the process 

9z . 1 d 2 z 2 \ dz 
dx 2 dx 2 J dx 



Az = ( — A + o ° 2 ) At + - ° Aw - ( 62 ) 



We can now obtain the Ito lemma from the path integral formalism. 

Let us interpret the function, z(x), as a change of variables from the variable, x, to the variable, z. The transfor- 
mation rules for the coefficients of the Langevin equation can be obtained by 

- 1 < 63 > 

k ' = £ k - m 

where 







h(x, t) = A(x, t) - i a(x, t) d x a(x, r), (65) 
h'(z, t) = A'(z, t) - i a'(z, r) d 2 a'(z, r). (66) 

The Eq. ([33]) gives directly the standard deviation for the new variable, z. Moreover, by using also the relation 

, dx d dz dx ( d 2 z dz \ 
dz dx dx dz \dx 2 dx x J ' 

we obtain 

Finally, if we compare this expression with 

*'=!*=£(A-^), (69) 

we obtain the drift for the variable, 2, 

., dz A 1 <9 2 z 2 
The expressions (^|) and ( f70| ) coincide with the results of the Ito lemma. 

IV. EXPECTATION VALUES 

The expectation value of a functional t/[x(r);r] on the stochastic process defined by the Langevin equation (|3^), 
with fixed initial and final conditions, x(t) = y, and x(T) = x, is formally given by 

< x,T I gr[x(r);r] | y,t >z = 

2?[VGx(t)] 5 [x(t);t] expj- / L[x(t),x(t);t] rfrl , (71) 
x(t)=y [ J 

The discretized expression corresponding to ( fn} ) is 

-00 /.+00 W+l 



• • • / n^? n (2^At)-»/ 2 Jote-utj-!) 
-°° »=i j=i 

x ff(x t ,xi, . . . ,x T ;i,*i, . . . ,T) exp{-5(xj,^;Xj-i,tj-i)}, (72) 



where ^(Xj, tj; Xj_i, tj-i) represents an 

appropriate discretization rule, while the function g(x t ,xi, . . . , Xy;t, ti, . . . , T) is a discretization of the functional 
<7[x(t);t]. In general, different discretizations could give different results. 

Obviously, the expectation value (|7l] ) can be also written in covariant form by substituting the Lagrangian, 
L[x(t), x(t); t], with the invariant one, £[x(r), x(r); t], by using the covariant discretization rule (^), and provided 
that also the functional, g[x(r); r], is discretized in a covariant form. 

A case of particular importance is when the functional can be written as 



£f[x(r);r] = exp V[x 



(r);r}drj. (73) 

Since the term, V[x(t);t], is invariant for a general transformation of variables, it can be simply included into the 
Lagrangian both in the non-covariant and the covariant formulations. Moreover, since such a term does not depend 
on the derivative of the stochastic variable, the differences among different discretizations of the functional (|7^) are 
of order larger than O(At), and, therefore, it can be discretized with any rule. 



10 



A. Mean value 



There are several examples of fvmctionals. The simplest one is the expectation value of the stochastic variable x at 
a given time t\. The functional can be formally written as 



(74) 



g[x(r); r] = / - t)x{t)<It t<t x <T, 

and the expectation value is given by 

<x,T\ x(h) \y,t> L = I dz p(x, T | z, ti) z p(z, ti \ y, t). 



For instance, in the simple case given by Eq. (|llj) the expectation value is easy to evaluate, and we obtain 

1 1 



<x,T\ x(ti) \y,t>i 



V^ttct 2 (T — ti) v / 2ncr 2 (ti - t) 
- +oc r (x-z) 2 



dz exp 



2 a 2 (T-tj.) 



z exp 



(y-z) 2 



2 a 2 (h - t) 



In absence of the integrand factor, z, we get the transition amplitude from t to T, 



p{x,T | y,t) 



1 



y/2ira 2 (T - i) 



exp 



(x - y) 



2 1 



2a 2 (T - t) 



(75) 



(76) 



(77) 



If we perform the integral with the integrand factor, z, we obtain the transition probability (|77j ) multiplied by the 
mean value of the functional, 



< x,T | x( tl ) | y,t > L = X{h h) p(x,T\ y,t). 



(78) 



B. Functionals for financial quantities 

Many financial quantities are defined as expectation value of functionals on a stochastic process, with fixed initial 
condition, x(t) = y. We will denote this kind of expectation value by the symbol, £x, which is defined by 



E L 5 [x(r);r] | x(f) = y =/ d n x < x, T \ 5 [x(r); r] | y, t> L 



(79) 



Note that this quantity could be formally written as a functional integral with only one extreme fixed, x(t) = y. In 
this case the functional measure would be completely invariant, and the whole functional integral could be formally 
written in full invariant form. However, from a computational point of view, this does not essentially simplify the 
procedure of calculation. 



1. Zero-coupon bond 
A first example is the evaluation of the quantity 



P(r t ,t,T)=E L 



f t T ,(r)dr | r(t) 



(80) 



which is the price at time t of a zero-coupon discount bond with principal 1, and maturing at time T . The variable, 
r(r), is the interest rate, and satisfies a SDE corresponding to the Lagrangian, L[r{r),r(T)\T\. Let us define 
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G(r T , T\r u t) = <r T ,T\e- C r(r)dT \r u t> L 

r r(T)=r T ( pT _ 1 

V[cr(r 7 r) r(r)] exp < — / L[r(r), r(r);r] dr 



where 



The price is then given by 



r(t)= 



L[r(r), r(r); r] = L[r(r), r(r); r] + r(r). 



/+oo 
dr T G(tt, T I ft, t). 
-oo 



(81) 



(82) 



(83) 



2. Capiat 

The price of a caplet is defined by the following expectation value: 



C{r u t,T,s) = E L 



r XT rMdr ( X - P(r T , T, *)) fl( x - P(r T , T, a)) \ r(t) = r t 



(84) 



which can be seen as the price at time t of a European put option maturing at time T on a zero-coupon discount 
bond with principal 1, and maturing at time s, where \ is the strike price. This expression can be put in the form 



C(r t ,t,T,s) = 



+ co p r.r(T)=r T 

dr T P[a(r,r)-V(r)] 

oo J J r(t)—rt 



( X -P(r T ,T,s))e(x-P(rT,T,s)) expj-^ {L[r(T),r( T );T]} dr 
dr T (x- P(r T ,T,s)) 0(x - P(r T ,T,sf) G{r T ,T \r t ,t), 



(85) 



where G(rx,T \ r t ,t) and L[r(r),f(r)]T] have been defined in Eqs. ( |Sl| ) and (g2J), respectively. This result has been 
reported, and explicitly calculated for the Vasicek model, by Jamshidian (cf. Eq. (11) in Ref. |^3| ). 



V. ANALYTICAL RESULTS 



In general, the explicit analytical calculation of a path integral is a formidable task, and it is possible in a very few 
cases only. A class of systems which allows exact path integration is characterized by quadratic Lagrangians [[14| [16| . 
In the next sections, we will consider some cases which belong or can be reduced to this class. 



A. Elementary cases 



1. Harmonic Lagrangian 



The harmonic Lagrangian is 



L[x,x;t] = x 2 + 



2^ 



(86) 



Note that this Lagrangian does not contain any term with a coupling between x and x, and cannot correspond to 
any stochastic process, but it is the starting point for the next calculations. The solution of the path integral for this 
process is a well known result, and it can be cast in the form 
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Iharmonic(x, T \ y,t) = e^-W exp{- ^g^ } , (87) 



where 



; (1 - e -^{T-t)\ 

° = " \l 2. ^ 



2. Harmonic Langevin equation 

Let us consider the stochastic process defined by the following SDE, 

Ax — —LuxAt + a Aw. 

The corresponding Lagrangian is 



(89) 



L[x, x; t] 



2 a 2 



\x + u> x] 



(90) 



This Lagrangian is equal to the previous one plus a coupling between the stochastic variable, x, and its derivative, x, 
and plus a constant factor. By using the prescription given in section IE, we get the following finite time transition 
probability, 



p(x,T \ y,t) = e 



x{T)=x ( ,T 

V[<j~ 1 x(t)} exp^ - / 

x{t)=y I Jt 



The path integral in Eq. ( |9l| ) is just Ihamonic( x i T I ?/>£)• Then, we obtain 



p(x,T | y,t) 



%/2tto- 2 



exp 



{ye 



j(T-t) 



xf 



2a 2 



where a is given by Eq. 



2a 2 2a 2 



dr 



(91) 



(92) 



B. The Vasicek model (Ornstein-Uhlenbeck) 

This model is very popular in the financial literature and it is defined by the stochastic equation |l7[ 

Ar = a(b - r)At + aAw, (93) 

where the variable, r, is the short term interest rate. The finite time transition probability, which, following the 
notation of Ref. will be denoted here by K(jt, T \ r t , i), is given by 



K(r T ,T\r t ,t) 



r(T)=r T ( ,T 

V[cT- x r{T)] exp ^ - / ' x 

r(t)=r t I Jt 



L v [r(T),r(r);T} dr\ , 



where 



L v [r,r;r} = [r-a{b-r)f - |. 



(94) 



(95) 



This case is similar to that of section V A 2 , with to = a, and x — > r — b. Then the path integral can be simply worked 
out by changing the variables (here we do not need the covariant formulation since this is a linear transformation of 
variables) . 



13 



1. Expectation value 

Actually, we are not really interested in the Green function for Eq. (|93|), but in the expectation value (|8l]), i.e. 



G(r T ,T\r u t) 



r(T)=r T 



r(i)=r t 



T>\a 1 r(r)] exp ■ 



L y [r(T),r(r);r}dT} , 



(96) 



where I/ V [r, r; r] = L v [r, f; r] + r. We can make the following linear change of variable, 



2 = r - b H — 5-. 

The Jacobian of the transformation is equal to 1, and the Lagrangian becomes 

fVr -i If , l2 a , u 2 i 
L [z , z; r}^—[ Z + az)} -- + &-_--. 



(97) 



(98) 



The last three terms can be integrated, and they give rise to a phase factor (we recall that we are using the mid-point 
prescription), then the integral (96) can be written as 



G(r T ,T \ r t ,t) = e 



-Ad 



z(T)=r T -b+a 2 /a 2 



z(t)=r t -b+(T 2 /a 2 

where the phase factor, Ad, is given by 



T>\a 1 z(t)] exp ■ 



1 ,■ \i a 
2^(* + «) "J 



dr 



A9 = 



r T - r t 



+ (T-t) 



b- 



2a 2 



(99) 



(100) 



Finally, since the remaining path integral is equal to that corresponding to the harmonic Langevin equation, by using 
the result (B2|) we obtain 



G(r T ,T\r u t) 



x exp 



V27TCT 2 

[{r t -b + a 2 1 a 2 ) e -< T -^ - (r T - b + cr 2 /a 2 )} 2 
2~d 2 



(101) 



with 



(1 - e -2a(T-t)) 

2a ' 



(102) 



This expression is the same that we find in the literature (cf. Eq. (13) in Rcf. |p~3[] ) . 



2. Zero-coupon bond and caplet price in the Vasicek model 

The quantity G(rx, T \ rt, t) allows to evaluate the price, P(r t , t, T), at time t of a zero-coupon discount bond with 
principal 1, and maturing at time T, by employing the formula (|83|). Since tt appears in a quadratic form to the 
exponent, the integral can be performed by completing the square and using the rules for Gaussian integration. The 
final result is given by 



P(r,t,T) = exp 



-B(T -t)r+ [B(T -t)-(T- t)](b - 



2a 2 



)-B{T-tf 



4a 



(103) 



where 



B(r) 



(104) 



A second task is to evaluate the caplet price given by Eq. (|85|). Such an expression has the structure of a truncated 
lognormal integration, and can be done analytically by usual methods. Here we do not discuss the details of the 
calculation which arc given in the literature. 
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C. The Black-Scholes model 



Let us consider the Black-Scholes model 



AS = a S At + a S Aw. 



(105) 



where S is the stock price. Here the metric depends on the stochastic variable, S, and we need a non-linear change 
of variables to perform the calculation. Therefore we must use the covariant formulation. 

The invariant Lagrangian corresponding to Eq. (105) can be obtained by Eq. (p>9|), with <r(S) = er S, and 

h(S) = a S — S. Since this is a one-dimensional problem, we can use the change of variables of Eq. fl58| ) which 

gives rise to a constant metric. Then the relation between the two coordinate systems is 



s = I logS. 



(106) 



Since the Lagrangian is invariant, we only need to exchange the quantities in the old coordinate system with the new 
ones, given by 



The invariant Lagrangian becomes 



Furthermore, the new measure is 



<j'(S')=<T, 

k'(S>) = §H { S) = (a-^). 

C[S',S';r} = ^[S'-(a-^)} 2 . 
V[<j- 1 S\t)] = S T P[c7(S')~ 1 S(t)]. 



(107) 
(108) 

(109) 
(110) 



Since in the new coordinate system the invariant Lagrangian (109) coincides with the mid-point one (see, section 
III A ), we can proceed in the usual way, and we get for the finite time transition probability, 



p(S T ,T\S t ,t) 



1 

St 



S'(T)=log(S r ) 



S'(t)=log(S t ) 



X>[<t _1 S'(t)] ex P { - 



2a 2 b + 2a 2 I 2 



dr 



where 



Ac/) 



logCSr) 2 2 f 

( a -—)dS' = (a-—) log(S t ) - log(5 T ) 

log(S t ) 1 1 L 



(111) 



(112) 



Finally, by using the formulae obtained previously, and collecting these expressions, we get the result of Black and 
Scholes. 



D. The Cox Ingersoll Ross model 

Another case where we need the covariant formulation of path integrals is the Cox Ingersoll Ross (CIR) model. 
This model has been developed to overcome a limit of the Vasicek model, i.e. the fact that the short term interest 
rate can take negative values. Obviously, this fact has not any real meaning. To avoid this problem Cox, Ingersoll, 
and Ross [[l8| have proposed the following process, 

Ar = a(b- r) At + a y/r Aw (113) 

(for a discussion of this equation, see, for example, Ref. [|l9|). Here the variance is proportional to r, and goes to zero 
with the short term interest rate. The choice of this particular dependence was related to the solubility of the model. 
As we will see, in this case the path integral ( |Sl| ) can be expressed in a closed form through an analytical function. 
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Let us start from the path integral written in covariant form. The invariant Lagrangian is given by the expression 

„2 



(p9|), with u{r) = a y/r, and h(r) — a (b — r) 

cr'(z) = a, 



If we introduce the variable z = we obtain 



h'(z) = % h(r) 



a{b - z 2 /4) 



a 

Jz 



Therefore, in the new variable, the invariant Lagrangian becomes 



C[z,z;t] 



1 
2^ 



a(b - z 2 /4) 



a 
4~z 



1 d_ 

2 ~d~z 



h'{z). 



(114) 
(115) 

(116) 



In order to evaluate the expectation value (81), we must add to the Lagrangian above the term, r = z /A, Furthermore, 
the measure in the new variable is given by 



If we now observe that the term proportional to z gives rise to the following phase factor (see, section IE), 



Acf> = r h'(z) dz = --Az T - z 2 ) + (2ab - log ^, 



with zt — and Zt = then the expectation value (|8lj) can be written as 



G(z T , T | z t ,t) = — exp I ^ 
zt " 



z(T)=z T 



:(t)=z t 



V[a 1 z(t)] exp • 



L e s[z, z; t] dr 



where 



L c s[z 7 z; t] = — - z 



2ct 2 



Cl 



2 , C -3 
C2Z H r, 

Z Z 



and the parameters, ci, C2, and C3, are given by 



ci = 



(7- 
,2 



C2 



C3 



1 

8^2 + 4' 

(4ab-3a 2 ) {Aab - a 2 ) 
8a 2 



(117) 
(118) 

(119) 

(120) 

(121) 
(122) 
(123) 



The effective Lagrangian, L ff, is not quadratic, but the term ~ — ^ is the same as the centrifugal contribution in the 

case of a two-dimensional harmonic Lagrangian, written in polar coordinates. Such a case can be solved in arbitrary 
dimension. The explicit expression and the algebraic and analytical manipulations of these rather cumbersome 
mathematical aspects are given in Appe ndix |^. 

The solution of the path integral (119) can be used to calculate the zero-coupon price by Eq. d83|). The final result, 
quoted in the literature (see, for example, fl9|| ), can be written as 



P(r t , t, T) — A(t, T) e - B (t> T )r t] 



(124) 



where the functions B(t,T), and A(t,T) are given in the Eqs. (C9), and (CIO). We stress further that in this case 
the anomalous 1/z 2 dependence masks the underlying quadratic problem. Once this connection is understood the 
solution is straightforward. 

The expression above represents a closed-form solution for the zero-coupon price in the CIR model. Obviously, the 
explicit knowledge of the Green function allows the evaluation of more complex functionals such as, for example, cap 
or floor prices. 

Moreover, the method can be applied to coupled SDEs. The general theorems on quadratic Lagrangians allow to 
solve by diagonalization a multi-dimensional coupled problem, the only constraint is the absence of higher order terms 
(0(x 3 ) or larger). 

Finally, this treatment can be extended to more complex problems with non-quadratic Lagrangians. The scientific 
literature suggests a variety of approximation techniques: for example, the perturbative, and the saddle point methods. 

In conclusion, the path integral approach, besides other analytical methods (see, for example, ]2p[|), represents a 
powerful tool of analysis. 
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VI. CONCLUSIONS 



We have described the path integral method as an alternative approach to find the solution of stochastic equations. 
We have shown that the method can be used in general cases, and that it is equivalent to the formulation in terms 
of SDE and PDE equations. The method is suitable in particular for defining functional mean values and path 
dependent problems. However, an exact analytical treatment is possible only in a few cases. A constraint necessary to 
get analytical solutions is the quadratic form of the Lagrangian describing the stochastic process. We have discussed 
in general this case, and the problem of the coupling between the stochastic variable and its derivative. Moreover, we 
have given the explicit solution for some important one-dimensional problems, and briefly discussed the generalization 
to quadratic multi-dimensional cases. Finally, we point out that the path integral treatment can be extended to non- 
quadratic cases by approximate analytical techniques, or by numerical methods. The numerical approach will be 
described in a forthcoming paper (see also Ref. pl[|). 



APPENDIX A: EQUIVALENCE BETWEEN THE PATH INTEGRAL APPROACH AND THE PARTIAL 

DERIVATIVE EQUATIONS 

We want to show that the short-time transition probability (|l3|) gives rise to the Fokker-Planck equation ( fTo| ) 
corresponding to the differential stochastic equation (ph. Let us write the following identity 



/-t-oo 
dz p(x, T + At | z, T) p(z, T \ y, t). 
-oo 

If we use the approximate expression (13), and employing the identity z = x + (z — x) = x + 77, we obtain 



/+00 
-00 



drj 



1 



exp 



(- 1] - A(x + t],T) At) 2 
2 At a 2 



p(x + r),T I y,t), 



(Al) 



(A2) 



where, for the sake of simplicity, we have taken a = constant. If we recall that r\ ~ 0(y At), and expand the above 
expression in 77 up to 0{At), we get 



p{x,T + At\y,t)~ I 



y/2-KAt er 



exp 



if 



2 At a 2 



x n-^A(x,T) 



r] 2 dA(x,T) x] 



(7 

d_ 

dx ' 



dx 



rj 2 d 2 



1 - 



At A(x,t) 2 



2a 2 



x \ p(x, T I y, t) + 77 — p(x, T\y,t) + Y ^ P( x > T \v^) 



Then, by performing all calculations, we find 

p(x,T + At\y,t)~p(x,T\y,t)-At 



dA(x, T) p(x, T I j^i) 1 9 2 a 2 p(x, T\y,t) 

2 



dx 



dx 2 



Since the L.H.S. of this equation can be also written, up to O(At), as 

p(x, T + At I y, t) ~ p(x, T I y, t) + At A p{x, T\y,t), 



(A3) 



(A4) 



(A5) 



by collecting all pieces together, we obtain the Fokker-Planck equation. 

The general case where there is an explicit dependence of the function, a, on x, and T, is more complex. We must 
expand the equation up to 0(r) 8 ), and 0(At 4 ). However after several cumbersome calculation we obtain the final 
result 



dp(x,T\y,t) _ dA(x,t)p(x,T\y,t) 1 d 2 a 2 (x, t) p(x, T \ y, t) 

2 Ox 2 



dT dx 
which coincides with the general Fokker-Planck equation ([[0] 



(A6) 
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APPENDIX B: EXPONENT EXPANSION 



In this appendix we describe a procedure to obtain an expression for the short-time transition probability which 
is correct up to any given order E^] . Here, for the sake of simplicity, we consider the one-dimensional case, but the 
formalism can be easily extended to arbitrary dimensions. Moreover we take a = constant, and A(x, r) = a(x). 

The starting point is that the transition probability must satisfy the Fokker-Planck equation (fiX\). Therefore, if we 
make the following ansatz, 



p(x,t + At | y,t) 



1 



\J2-k Ata 2 

the function, f(x,y, At), must satisfy the equation 
d 



exp 



~(x-y) 2 

2 At a 2 



f(x,y,At) 



(Bl) 



dAt 



f{x, y, At) = -a{x) f(x, y, At) + \^^ /0, V, At) 

_ A a (x) + ( x -y^ JL f( x v At) ~ {X ~ V) a(x) 
dx [)+ At dx n ' V ' ' Ata 2 1 



Let us now expand f{x, y, At) in powers of At, 

f{x, y, At) = f (x, y) + At f x {x, y) + At 2 f 2 {x, y) + At 3 f 3 (x, y) + 
then, by substituting in (JB2|), we get the following set of recursive equations, 



(B2) 



(B3) 



d x fo{x,y) 
fi(x,y) 



— I a ( x ) 



= -(x - y) d x fi{x,y) + \ d x a(x) + 

2f2(x,y) = -(x-y) d x f 2 (x,y) + ±o 2 d 2 Ji{x,y) 

3/s (a, 2/) = -(x-y)d x f 3 (x,y) + \ a 2 d 2 f 2 (x, y) - 

4:h(x,y) = -(x - y) d x f 4 (x,y) + | a 2 d 2 f 3 (x,y) - 



a{xf 
2a 2 



(B4) 



1 „2 



^_a 2 {d x h{x,y)) 2 
he 2 d x fx(x,y) d x f 2 (x,y) 



The first equation is decoupled by the others and defines an overall phase factor. The term fo(x,y) is related to the 
integral of a(x), and it is just the phase factor (p8|). The second equation gives the first order approximation in At; its 
value can be used as input for the successive equation, and so on. Therefore, by a simple iteration, we can compute 
the function f(x,y,At) up to the required order. In general, the structure of these equations is 



nf n (x,y) =-{x-y) d x f n (x,y) + W n {x,y), 
where W n {x, y) depends on the functions, fi(x, y), . . . , f n —i(x, y). The solution is given by 



(B5) 



f n (x,y) = / d£ £ 
Jo 



Actually the integral appearing in Eq. 
further expansion in Aa; = x — y |E3| . 



" - 1 W n (y + ax-y),y). (B6) 

) can be rather complicated, but the problem could be simplified by a 



APPENDIX C: EXPLICIT CALCULATION FOR THE CIR MODEL 



The path integral (119) has an explicit solution given by (see, for example, ]16|) 

7 y/Z T F t 



G{z T ,T\z u t) 



2 JA01 
exp < > 



x exp 



7 
4ct 2 



z t I a i 2o 2 sinh 
(z 2 +z 2 )coth{l(T-t) 



:(T-t) 



z T z t 7 



2a 2 sinh 



Q(T-t)) 



^a 2 b(T-t)/<7 2 



(CI) 
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with 1^ (x) the modified Bessel function of index \i, 



Va 2 + 2a 2 , 



and 



! c 3 2ab 



1. 



(C2) 



(C3) 



In order to get the formula for the zero-coupon price this expression must be integrated over the final variable, tt'- 

p +00 



P(r u t,T)= dr T G(z T ,T\z t ,t)= dz T -j- G{z T ,T \ z t ,t) 
Jo J z 



The integral looks rather formidable, however it can be handled thanks to the integral 



+00 



e -a*V +1 Iptfz) dz = 



(2«) 
In the equation above the parameter a, and (3 are given by 



exp 



4a 



and 



7 

a = — 77 coth 
4cr 2 



/3 : 



7 



(T-t) 



a 
4^2' 








2cr 2 sinh 





Therefore, by performing the integration, and inserting the original variable, r t , we obtain the expression 



P(r t ,t,T) = exp 



4q: 



f ar * 7H [7, T 



7 



A a 2 a sinh 



exp 



-(T-t) 



Note that the dependence of P(r t ,t,T) on r t is only due to the exponential, exp[— B(t,T) r t ], where 



B(t,T) = 



7 



4aer 4 sinh 2 
Finally, if we define A(t, T) by 

A(t,T) 



liT-t) 



a 7 
5 + 4r coth 



2(T-t) 



2( e 7(T-t) _ !) 



(7 + o)(eT( T -*) -l)+2 7 ' 



7 



4 tr 2 a sinh 



r(r-i) 



exp 



2 7 e 



7 (T-t)/2 



(7 + a)(e7( T -*)-l) + 2 7 



2a6/o 



,a 2 b(T-t)/cr 2 



(C4) 



(C5) 



(C6) 



(C7) 



(C8) 



(C9) 



(CIO) 



we can cast the expression for the zero-coupon price, P(r t , t, T), in the form given in Eq. (124). 
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